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We review recent progress toward understanding of sQGP. The phenomenological part 
includes discussion of elliptic and conical flows at RHIC. Then we proceed to first quantum 
mechanical studies of manybody states at T > Tc, the "polymeric chains" q.g.g...q and 
baryons. A new model for sQGP is a classical dynamical system, in which color vector is 
changed via the Wong equation. First Molecular Dynamics (MD) results for its diffusion 
and viscosity are reported. Finally we speculate how strong correlations in matter may 
help solve puzzles related to jet quenching, both the magnitude and angular distribution. 

1. Introduction 

The idea that QGP at RHIC, T = (1 — 2)T C , seems to be in a strongly coupled regime 
(sQGP) was introduced in 2003 [d 121 El and was discussed a lot at the "discovery" BNL 
workshop [Ej and QM04. Let me start with the simplest physics point possible. Although 
at T > T c quarks are deconfined, the energy needed for separation of quarks close to T c 
is huge (see FigGK up to U ~ 4GeV. One implication is simple: the charges simply 
cannot get separated until very high T. The second: as the ratio U/T ~ 10 goes into the 
Boltzmann exponents, any perturbative approach is completely hopeless. 

Instead of going to many other arguments, let me simply list main reasons which de- 
mand such a radical change, from traditional weak coupling to strong coupling methods: 

1. Collective phenomena observed at RHIC lead the hydro practitioners to view QGP as 
a "near perfect liquid" [ [TJ E] ; 

2. Feshbach-type resonances due to marginal states may lead to large cross sections 

3. Classical e/m plasma can be a good liquid too, if sufficiently strongly coupled. 

4. A close relative of QCD, the J\f—4 supersymmetric Yang-Mills gauge theory can be 
studied in a strongly coupled g 2 N c — > oo regime at finite T via the AdS/CFT correspon- 
dence: the results are very close to what we observe at RHIC. 

Since I cannot cover all these topics in depth, let me just provide a more detailed list 
of these arguments, with at least some statements and references: 

(la) The departure of elliptic flow data from the hydro prediction happens only at rather 
high p t ~ 1.5 — 2GeV, from which the estimated viscosity-to-entropy ratio rj/s = .1 — .2 

is more than order of magnitude lower than in pQCD. 
(lb) Another transport coefficient, the charm diffusion constant D c deduced from single 
electron Raa and t>2 (much discussed at this meeting), is also an order of magnitude lower 
than pQCD estimates [[7j. 
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Figure 1. (a) The separation energy U(T,r — > oo) — U(T,r = .2fm) in GeV vs T/T c , 
calculated from the free energy [ |3] by removing the entropy term, (b) Dependence of 
various states' binding energy (in units of T c ) on the temperature. 



(lc) New hydrodynamical phenomenon suggested recently the conical flow, maybe 
explains why secondaries from a quenched jet fly preferentially to large angle ~ 70 degrees 
consistent with the Mach angle for (time-averaged) speed of sound. 
(2a) Combining lattice data on quasiparticle masses and interparticle potentials, one in- 
deed finds a lot of bound states [E] and resonances [[TUj at T > Tc 
(2b) The same approach explains why n c , J/ip remains bound till near 3T C , as was directly 
observed on the lattice [HI]; 

(2c) this approach was experimentally demonstrated to work for ultra-cold trapped 
fermionic atoms Li 6 , turning it to near prefect liquid as well when the Feshbach reso- 
nance leads to the scattering length a — > oo. Experiments on oscillations found a sharp 
minimum in damping near the resonance, reducing it by about 2 orders of magnitude. 
According to our study [E2] two lowest modes are well described by hydro with "quantum 
viscosity" rj ~ 0.37m, in a way nearly as low as that of sQGP. 

(2d) Heavy-light resonances in sQGP can explain the value of the charm diffusion con- 
stant [ITS]. 

(3a) The interaction parameter T =< potential energy > /T in sQGP is obviously not 
small, O(10). At such Y the classical strongly coupled e/m plasma is a good liquid: we 
find the same in our classical version of sQGP as well; 

(4a) The EoS for finite T J\f=4 SUSY YM is similar to what is seen on the lattice in the 
RHIC domain, namely p/p ideal = [(3/4) + 0((g 2 N c )- 3 / 2 )} [UJ; 

(4b) At infinite coupling there is a finite limit of viscosity n/s => 1/47T again close 
to the RHIC value. The "R-charge" diffusion constant is low as well. 
(4c) One may even think about AdS/CFT "gravity dual" to the whole RHIC collisions, 
with black hole production [^] which then fly away from the test brane [I17j. 
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2. Systematics of collective flows 
2.1. Elliptic flows 

Collective flows observed at SPS and RHIC are quite accurately reproduced by ideal 
hydrodynamics for the QGP/mixed phase, complemented by a hadronic cascade for the 
hadronic phase [EI]. The radial flow predictions are in quantitative agreement with the 
data for all secondaries, from n to fi~. Elliptic flow is described by v% =< cos(20) > 
depending on 6 variables V2(s,pt, Mi,y,b, A), the energy, transverse momentum, particle 
mass, rapidity, centrality and system size. Hydro works for not-too-large p t < 2GeV 
(above it a different and still poorly understood regime starts.) 

(i) The energy dependence was the issue I discussed already at QM99[HH], v 2 (s) was 
predicted to be smoothly rising from SPS to RHIC by about factor 2, as stiffer EoS of 
the QGP replaces that of hadronic gas and mixed phase, see more in [EI]- It was soon 
confirmed by the first RHIC data 1 . 

(ii) the (pseudo)rapidity dependence v 2 (r/) has a triangular shape (PHOBOS). Although 
the first 3-d hydro by Hirano et al were not able to reproduce it, now with a hydro + cascade 
approach [|20] it works well. The physics is basically the same as for V2(s): in fact v 2 dis- 
plays good "limiting fragmentation", it depends mostly on local density v 2 (dN/dy(s,y)) 
rather than on s and y by itself. 

(iii) The u 2 (pt,Mj) dependence for different species, known as the "fine structure" of 
the elliptic flow, was discussed here by R.Lacey [ I2T] who pointed out a scaling relation 
v 2 ~ Miyl where y t is the transverse rapidity. Suggested by hydro, it is perfectly satisfied 
by the data for all secondaries. 

(iv) Centrality and size dependence: ideal hydro is scale invariant and v 2 is basically 
given by the spatial deformation, thus v 2 (b, A)/e 2 (b, A) « const (b, A). This relation is in 
excellent agreement with the data. 

Remaining hydro-skeptics include Bhalerao et al [ 22], who ascribe the rise in v 2 (s) (or 
that in v 2 (rj) toward midrapidity) to an "incomplete equilibration". It was pointed out 
earlier by H.Heselberg et al that hydro scaling (iv) should be violated for sufficiently 
peripheral collisions, by a "dilute regime" . Unfortunately, due to experimental difficulties 
that regime was never really clearly observed. Bhalerao et al pointed out possible depen- 
dence on the system size, predicting violation of hydro scaling v 2 (CuCu) / v 2 (AuAu) ~ 1/3. 
Recent CuCu data do not agree with them while the hydro scaling is still satisfied. 

An overall blue sky still has two clouds: (i) v^jv\ seem to be a constant, but not .5 
as ideal hydro predicts, but about 1.2. It may be related to the issue of v 2 fluctuations 
discussed by Mrowczynski and myself [121].). (ii) the HBT radii remain a puzzle. Para- 
doxically, we can better describe early stages (when the elliptic flow is formed) than very 
late dilute ones. It may well be that our freezeout conditions via naive cascades are no 
good: effective potentials, resonances etc may be different than we thought 2 . 



^^Hydro with fixed freezeout, at Tf, lead to a non-monotonous curve which is not the case. 
2 One recent attempt, by Cramer and G.Miller, has generated discussion here. I complained that they do 
not account for the pions absorbed by their optical potential on the way out. B.Muller, in a summary, 
defended them because "the source term in the HBT formalism describes the vertex position of the last 
inelastic interaction". As the author of that formalism (in 1973) I cannot agree more: thus Cramer and 
Miller better follow their pions to that last interaction. 
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2.2. Conical flow 

Let me start reminding a general kinematics of the most elementary process 1 => 2, 
possible for a jet in matter. No change in the mass implies = (p — k) 2 —p 2 = —2(pk) + k 2 . 
Provided the last small term (recoil) can be dropped, the angle between k and p is then 
cos{6) = (k p /kp), For very fast particle po/p = 1 and thus the condition for the process 
is simply k < k. 

Now, the RHIC data on 2-particle correlation with a triggered jet show that most 
radiation goes into a peak at angles 60-70 degrees, so that cos{6) « 1/3 and thus that 
must be the ko/k ratio for what is emitted. This fits very well to sound waves, for which 
this ratio is the speed of sound c s , and its average value over the duration of the collision 
is indeed about 3 1/3. One simple way to test that was suggested by Antinori and myself 
[I25j: if the jet is a b quark (which can be tagged experimentally) Po/p = 1/v > 1 which 
causes the cone to shrink, till it goes to zero at the critical velocity v — c s . Gluon radiation 
behaves in the opposite way with decreasing v, and never shrinks to zero. 

A number of people suggested a possibility of gluon 
Cherenkov radiation, see the talk by Ruppert for de- 
tails. There is very little support in data or theory 
for that: but even the optimists find the angles for 
Cherenkov emission of gluons with p ~ 1 GeV to be 
at least order of magnitude smaller than the observed 
peak at ~ 70 degrees. 

Whether the mechanism of dE jdx is a gluon radi- 
ation or scattering or ionization, a small mean free 
path implies rapid local deposition of energy, with 
subsequent hydrodynamical process. This issue is 
discussed in detail in the talk of my collaborator 
Casalderrey-Solana [|H]: The observed shape of the 
peak, its magnitude and p t dependence is in rough 
agreement with the data. The issue of flow subtrac- 
tion seem to be resolved (see B.Cole's talk), but 3- 
particle correlations are still too controversial to com- 
ment. Let me only remind that our hydro solutions 
are obtained with arbitrary viscosity, and that the observed effects already constrain the 
viscosity more than all data on v% since the gradients are larger for conical than for elliptic 
flow. 

3. First quantum manybody states above T c 

The existence of bound states of quark and gluon quasiparticles at T > T c , suggested 
in [|2] was much discussed at this meeting. Let me first report on recent studies [ 26J, 
done in collaboration with grad. student J.F.Liao. 

One result is about polymeric chains of the type qgg.-gq which appear naturally in a 
string picture 4 . For this problem the most useful coordinates in this case are not the usual 




Figure 2: Schematic picture of 
matter, as the temperature grows 



3 Note, we don't mean here speed of sound squared! 

4 Recall that a gluon has two color indices and can thus be connected to two strings. 
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Jacobi but "chain coordinates" , using which we have proven that such polymers have the 
same binding energy per bond? as for qq mesons. We also studied (variationally) two 
3-body problems: the closed (3-) chains of gluons (ggg) and (qqq) baryons. Their binding 
energy (and range) is shown in Fig[j"] one can see that the former one is quite robust 
while baryons have quite marginal binding (because of smaller relative charge of quarks). 
Before we go forward with a general discussion, let us try to summarize the proposed 
scenario as a single picture, see Fig^b). From relatively short string-like configuration 
of color fields at low T, fig (a), one moves to longer strings (b) at the critical point [125]. 
New is picture (c) which depicts "polymeric chains" considered in this work, significant 
at T = (1 — 1.5)T C . Eventually, at high T, one goes into (d) with independent quark and 
gluon quasiparticles, neutralized by isotropic Debye clouds. 

Let me now briefly touch on criticism of the 
"bound state QGP" raised at this meeting. First, 
let us be clear what is being criticized. The bound 
states are not ad hoc ideas but the direct mathemati- 
cal consequence of: (i) quasiparticle masses from the 
lattice [I2H|; (h) lattice-based effective interaction be- 
tween colored charges [[""] ; and (iii) quantum mechan- 
ics. None of the critics ever suggested that either (i) 
or (ii) or (iii) may be wrong. Where the direct test 
has been made, for the charmonium states [E], one 
finds a good agreement. And it does not take much 
to realize that quasiparticles with M ~ 1 GeV should 
behave like the charmed quarks. 

V.Koch et al [I3U] and F.Karsch[ I3T] argued that 
they do not see a significant contributions of di- 
quarks (and baryons) in the lattice data on d n (T) = 
d n (p/T 4 )/d(fi/T) n \ fl=0 with n=2,4,6 at T > 1.2Tc. 
It is hardly surprising, as there are not so many states 
of those and the expected masses are rather large, 
2M q and 3M q , respectively 6 . (In [[""] we have not 
even included diquarks and baryons in pressure, to 
which they clearly contribute too little.) Karsch [ 
I3T] further argued that the ratio d i /d 2 ~< B 2 > 
gives directly the mean constituent baryon number 
B, and since this ratio gets close to 1 above T c , all 
bound states with B > 1 are excluded. But the same 
reasoning gives d§/d± m< B 2 > as well: this ratio is 
however not close to 1 but ~ —10. Does it exclude 
any quark gas model as well? 

In fact the argument is simply too naive and ig- 
nores a lot of things. What is worse, Karsch et al 

5 Not per particle, so long enough chains have twice more binding per particle than mesons. 

6 In fact all their arguments can be repeated verbatim for say T ~ 50MeV or so, with the same conclusion: 

no visible baryons in a hadronic gas. They do exist as states there, of course. 
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have not explained the most prominent features of his 
data, the peak in d±(T) and the "wiggle' in d 6 (T). see Fig|21 or large flavor-nondiagonal 
d 4 {T),d 6 {T). 

It was pointed out in ref.[|33|, that since the quasiparticle masses may depend on fi 
their derivatives such as M" = (d 2 M/dfi 2 )\^ = o must be included in all formulae. The 
same is true for masses of the bound states such as baryons. It is shown in our paper 
[ that all above mentioned features of higher susceptibilities are naturally explained 
by the expected change in the baryon (N,A) mass M B (T,fi). At the phase boundary 
it is expected to grow from the "vacuum mass" ~ lGeV to much larger value in sQGP, 
rs 3M q . As a result M'g is large and changes sign at the inflection point near T c : that is 
why there is a peak in cLt(T) and the "wiggle' in d 6 (T), see the baryon curve b in Fig 01 
(The effect of binding of qq and qg states is not included in this plot: in the latter case 
it can increase its contribution by factor 2-3 and make a better agreement with data, 
especially for d 2 .) 



4. Strongly coupled colored classical plasma, studied via Molecular Dynamics 

Nice introduction to strongly coupled Abelian e/m plasmas has been provided here by 
M.Thoma, so let me jump directly to the point. The interaction parameter is defined as 

r = C c a s n l ^N corr /T (1) 

where the color Casimir is C c = 4/3 and 3 for q and g. The quasiparticle density 7 is 
proportional to large number of degrees of freedom in QGP 7iV3 ~ N^T where N dof = 
16 + 2*2*3* Nf ~ 50. The number of "strongly correlated partners" N corr is just 1 
for binary states, 2 for polymeric chains we discussed above, but goes up to 12 for a fee 
crystal (see below). Combining all, even for a s ~ 1 one finds T > 1, possibly up to O(10). 
Further, the structure is exponentially sensitive to T as it goes to Boltzmann factor. And, 
as we show below, T ~ 10 is precisely where the plasma liquid is most perfect. 

New model B.Gelman, I.Zahed and myself are proposing [133] is based on the following 
main assumptions: (i) The particles are heavy enough to move non-relativistically, with 
masses 8 M » T; (ii) Their interaction is dominated by colored electric (Coulomb) inter- 
actions, with all magnetic effects (spin forces, etc) ignored; (iii) Their color representations 
are large, so that color operators t a can be represented by classical vectors. 

Dynamics of color vectors as well as particle coordinates are described by classical 
equations of motion (EoM) 9 known as Wong equation 10 . 

dC a 

m ^ = 9f abc P^C c (2) 

7 Quasiparticles inside various bound states should also be included. 
8 Rccall that close to T c the M/T ratio for quasiparticles reaches 4-5. 

9 Anothcr model, amenable to MD simulations, has been proposed at this meeting by the Budapest group 
[ lrifi|: but the color part is not classical and is subject to complicated stochastic transitions. Not only 
color is not conserved in their model, but (as far as i can tell), the energy cannot be conserved either. 
10 It can also be rewritten as a set of canonically conjugated equations, for one x-p pair for SU(2) and 3 
pairs for SU(3), see [|27| for refs. B.Muller in summary expressed doubts whether this model is gauge 
invariant: no reason to worry: it obviously is, this is how Wong derived it. 
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where is a local field induced by other particles. The total color and total energy 
are thus conserved (which is monitored during actual MD calculations). The interaction 
potential is proportional to the dot product of the unit color vectors times the potential. 
The latter eventually should be deduced from lattice simulations, we so far use just a 
Coulomb. In order to stabilize the system, we have also added a short-range repulsion 
which is suppose to mimic the quantum-mechanical localization energy V s h ort = h 2 /Mr 2 , 
like it is sometimes done to describe systems like solid He. 

The main variable parameter of the model is obviously the temperature T. Eventually, 
in applications to sQGP all 4 parameters mentioned above are in fact some functions of 
T, so such applications of the model in the narrow sense would only be restricted to some 
1-d line in a 4-d parameter space. But in order to understand them better it would be 
necessary first to describe the properties of the model in a wider sense, in all its parameter 
space. 

The closest physical e/m analog of our model is the classical two-component ionic 
plasma, e.g. a molten ionic salts 11 such as NaCl. For large r — ► oo this system 
freezes into the fee cubic lattice with alternating positive and negative charges. For 
smaller 1 < V < T c ~ 80 one deals with a strongly coupled ionic liquid. In the non- 
Abelian case one also expects that the system gets frozen at very large T, with the same 
cubic crystal plus quasi-Abelian ferromagnetic (alternating) order of the color vectors. 

Molecular dynamics (MD) simulations can follow 
evolution of classical systems, weak or strong cou- 
pling. Our studies of that are at exploratory stage 
as of now, done so far with 4 3 and 6 3 particles with 
the SU(2) color vector (a vector on a 3-sphere). In 
the QM05 talk I have shown some movies of melting 
a crystal by "heating" a system with a small random 
force: a bit of friction can cool it back. Even when 
pictures look homogeneous and the color vector looks 
isotropic, there are local correlators of neighboring 
particles at large T. 

Our interest is not so far in detailed studies of its 
microscopic structure per se, but in collective ex- 
citations - the phonons and plasmons, as well as 
their widths reflecting transport properties such as 
the (self) diffusion coefficient D and viscosity r) by 
Kubo formulae 12 . With increased potential energy 
(or decreased T) the particle mobility decreases, to- 
gether with D. But viscosity has a more interesting 
trend, with a minimum and subsequent rise toward 
a "glassy" liquid and solid regimes. This happens because the momentum transfer can 

11 It is a good approximation to think an electron to be completely transferee! from Na to CI and T still 
low enough not to excite other electrons. 

12 In fact MD is well suited for such studies, while in lattice QCD the Euclidean time makes it next to 
impossible to use Kubo-like formulae and access such transport properties. 




Figure 4: Our first MD 

data for the SU(2) classical 
plasma: the particle diffusion co- 
efficient D and viscosity rj ver- 
sus T =< \Pot.energy\ > / < 
Kin.energy >. 
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be achieved not only by particle propagation but also by phonons/plasmons, which have 
large mean free path in perfect solids. Thus an optimum exists, with "the most perfect 
liquid": as FigEI shows, it is at T ~ 10 as well. 



5. New views on jet quenching and charm trapping 

The relevance of substructure of sQGP should be even more important for transport 
properties. It was first suggested by Zahed and myself in [|2| that existence of marginal 
states must increase rescattering and thus dramatically reduce viscosity (mean free paths), 
leading to a collisional (hydrodynamical) regime of expansion. Since different states get 
marginal at different T, one may hope this mechanism to work at all T up to about 2T C , 
the highest temperature at RHIC. 

We will now argue that the situation can be different for jet quenching. Its radiative 
theory tells that all properties of matter come one combination, q = 5 — 15 GeV 2 / fm, 
the values currently under discussion. In most papers this value is expressed as the gluon 
density, with dN/dy = 1000 — 3000, and some people argued this may be too large and 
contradict to entropy produced, see e.g. [EOj. However the relation between q and dN/dy 
assumed is simply invalid in the sQGP regime. (Think of a iVaC7-like local structure, with 
alternating charges: the electric field between the ions is a coherent sum of fields from 
positive and negative charges, which increases local fields and q compared to randomly 
placed charges.) 

Another way to see why jet quenching (and charm equilibration) should be enhanced 
near T c is to think about collective momentum/energy sharing in multibody correlated 
clusters we studied above. In particular, an admixture of polymeric chains is known in 
material science to be very effective mechanism of the momentum distribution over larger 
volume of matter. A famous example is Kevlar fibers added to epoxy (or other plastic), 
now applied in wide range of applications, from tires, boats etc. to such exotic ones as 
"bullet-proof vests" and even "anti-mine boots" . 

If this is so, there are consequences which can immediately be checked in experiment. 
If indeed only a "well correlated" (polymerized?) sQGP provides large dE/dx, it should 
only happen at T rather close to T c . Since at RHIC the initial temperature is about 2T C , 
jets may propagate with smaller losses, till the matter cools down and correlates properly. 
Such a delay would affect the angular distribution of the jet quenching phenomenon for 
non-central collisions. 

In fact the original idea of "jet tomography" via jet quenching is in serious trouble for 
quite a while, because the most natural assumption - the energy loss proportional to the 
matter density - is in strong contradiction with the observed strong angular dependence 
of jet quenching [EH], predicting too weak azimuthal asymmetry for non-central collisions. 
It was however recently pointed out by Pantuev [EH] that a better description of data can 
be achieved if the jet quenching at the highest RHIC energy is switched on after some 
"latent time" of about 2.2 fm. This time quite reasonably matches a cooling time from 
T ps 2T C toTsi 1.5T C at RHIC. 
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6. Status of the theory 

Studies of sQGP are developing very rapidly, with many new ideas and connections to 
other fields emerging daily. Instead of a summary, let me comment instead on the rather 
complex status of QGP theory in general. 

Before RHIC there were many approaches which looked reasonable, but then RHIC data 
put it under severe tests. Most of them have not survived, failing to explain collective 
flows, especially the elliptic t> 2 13 . It took the courage of Gyulassy and Molnar [EH] to 
tell us that a parton cascade can only get close to the data if the cross sections are huge, 
<jgg ~ 30mb. But if so, the number of particles interacting at any moment is large and 
the very concept of a cascade - free motion interrupted by scattering events - looses 
its meaning 14 . Lattice studies have found very large potentials between colored charges 
near and above T c . Now a challenge to all of us is to learn how such a strongly coupled 
system of quasiparticles can be organized. The paradigm shift, from wQGP to sQGP, had 
definitely occurred. Of course, all these events created a lot of stress in theory circles. 

Is the lesson learned? Well, a lot of people are in denial. Example: in the theory sum- 
mary B.Muller, after briefly acknowledging a discovery of "near perfect liquid", proceeds 
to "new ideas" such as a parton cascade with plasma instabilities. This subject, reviewed 
by Mrowczynski, is indeed interesting, and recently numerical studies have shown how 
such instability may develop in a quasi-Abelian regime. But now we are approaching 
its moment of truth: its time to do simulation without boxes, with appropriate initial 
distribution of partons, testing whether the instabilities can indeed lead to a "little bang" 
and i>2- I may be wrong, but there are many reasons to think that a collisionless weakly 
coupled plasma will fail this test. How large gluon masses and potentials seen on the 
lattice can be put into negatives by feeble magnetic effects? How can those explain heavy 
charm quark stopping/flow, much discussed at this meeting? 

The situation is well illustrated by a cartoon shown in the summary by B.Muller 15 , 
displaying a chicken afraid to cross the street. Yes, crossing from familiar wQGP to the 
sQGP is scary. My advice: don't be a chicken, learn to fly, be an eagle. High above there 
is no fear, and one can see what people do in other fields. And, last but not least, there 
is basically no alternative: one can sell non-working theories only for so long. 
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